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INTRODUCTION 

The following discussion is one of the results of the research for determining the vertical 
structure of the wind-driven circulation in an enclosed basin comparable in size to the 
Pacific Ocean. 

The first attempt to verify the effect of prevailing winds in maintaining the oceanic 
circulation was undertaken by Sverdrup (1947). According to his result, the oceanic cur¬ 
rents in the eastern part of the Equatorial Pacific are largely fed by the energy of the 
winds blowing over the surface of the ocean. Reid (1948) confirmed this conclusion. At 
nearly the same time, Stommel (1948) could explain the intensification of the wind- 
driven currents along the west coast of an ocean by assuming the existence of horizontal 
friction and the meridional variation of Coriolis force. Altogether, these investigations 
have enabled us to ascribe the major part of the oceanic circulation to the result of the 
superincumbent wind system prevailing over the oceans. 

Munk (1950) published a very important paper on the wind-driven circulation of the 
oceans. The next year Munk and Carrier (1951) treated the circulation of the North 
Pacific, regarding this ocean as a triangular basin. They could explain the pattern of the 
actual ocean circulation very well, so that there seems to be little left to discuss on this 
subject, as far as the major characteristics of the general circulation in the Pacific Ocean 
are concerned. 

Hidaka (1951) solved the problem of the general circulation which would be produced 
by both zonal and anticyclonic wind systems. In this computation, spherical co-ordinates 
were used, thus taking the effect of the sphericity of the earth into account. But both 
assumptions of the wind distribution gave no essential difference in the results except for 
the magnitude of mass transport. Moreover, the result for zonal distribution gave no 
sensible difference in the pattern of the circulation compared with Munk’s which was 
derived by using a rectangular co-ordinate system, except in the magnitude of the mass 
transport. These facts show us that we have only to treat the circulation driven by a 
zonal wind system. 

An earlier paper (Hidaka, 1950) contained a theory of ocean circulation using the cur¬ 
rent velocity in place of the mass transport. The analysis was complex because the vertical 
variation of the movement of water had to be taken into account. The result was, none¬ 
theless, quite ridiculous; no perceptible concentration of the streamlines toward the west 
coast could be found. The explanation for this result was that the approximation to the 

1 Contribution from the Department of Oceanography of the Agricultural and Mechanical College of 
Texas, Oceanographic Series, No. 42. Based in part upon work done under the sponsorship of the Office 
of Naval Research and the Bureau of Ships. Manuscript received July 2, 1953. 

2 Geophysical Institute, Tokyo University, Tokyo, Japan. Visiting Graduate Professor, 1952-53, Depart¬ 
ment of Oceanography, Texas Agricultural and Mechanical College. 
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y = -“j a (60° S) 

Fig. 1 . A rectangular ocean comparable in size to the Pacific Ocean. 


mass transport was inadequate for its east-west variation. Thus, the solution smoothed 
out the western currents and boundary vortices which were apparent in Munk’s paper. 

The mass transport method which has been adopted by Defant (1941), Stockmann 
(1945-46), Sverdrup (1947), Reid (1948), Munk (1950, 1951), and recently by Hansen 
(1952) is surely eminent, especially in that it enables us to reduce the analysis to two 
dimensions and makes the mathematical procedure very simple. Moreover, it is not 
necessary for us to consider the vertical variation of density and vertical coefficient of 
eddy viscosity. These authors have indeed contributed greatly to the solution of many 
important problems of oceanic circulation by this method. The author himself also em¬ 
ployed this method several times in discussing the problems in this direction. However, 
it is impossible for this method to show how the wind-driven circulation varies in a vertical 
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direction. Neumann also expressed recently (1952) his opinion as to the necessity of a 
dynamic treatment of ocean currents as a three-dimensional problem. 

All these circumstances lead one to recompute the general circulation of the Pacific 
Ocean under these modified conditions and assumptions. The present investigation is one 
of the results of the author’s efforts in this direction. We here treat the general circulation 
of the water in a square ocean comparable in size to the entire Pacific Ocean basin. 
Spherical co-ordinates, which were used in a preceding paper (Hidaka, 1951), are not 
used here, partly in order to avoid mathematical difficulties but mostly because the two 
systems of co-ordinates did not give any important difference between Munk’s and the 
author’s results except for the magnitude of mass transport. The value of the lateral 
mixing coefficient is taken as 3.08 X 10 7 c.g.s., a value consistent with the research of 
former investigators. The wind system is considered zonal, because this assumption is 
far simpler for the subsequent analysis, and also because no essential difference has been 
found between the results obtained under the assumptions of zonal and anticyclonic wind 
distributions. Of course, the variation of the Coriolis parameter with latitude is taken 
into account. The use of current velocity in place of mass transport makes the mathe¬ 
matical analysis many times more complicated because the problem is now three dimen¬ 
sional. But the result will be of importance because it should give an idea of the vertical 
structure of the wind-driven circulation of the oceans. 

THEORY 

The dynamic equations of the stationary ocean currents, taking both vertical and 
horizontal mixing into account and neglecting the nonlinear terms, are 

. d 2 u d ( . du\ dp 

- 4, V + SV 4 -*J + 2 " s,n *'”’-to = 0 ' (1) 

, d 2 v d ( . dv\ 0 dp „ 

A 'W‘ + di\ A -ai) ~ “ Ty = 0 

where u and v are components of the current velocity in x (eastward positive) and y 
(northward positive) directions, p is the pressure, p the density, Ai and A z the horizontal 
and vertical coefficients of eddy viscosity of the water, oj the angular velocity of the earth, 
and <f> the geographical latitude. The axis of z is taken positive downward, the origin 
being placed on the undisturbed sea level. 

The boundary conditions to be satisfied on the surface (z = 0) and at the bottom 
(z = h) are 

z = 0: —A z ~ = t x (x, y); —A z = t v (x, y) (2 

and 

z — h: u = v = 0. (3) 

Here both coefficients of eddy viscosity are supposed to be constants. The conditions to 
be satisfied along the coasts are also necessary. These are simply that there is no water 
flow across and along these coasts. If the coasts consist of vertical cliffs, we have 

u = v = 0 (4) 

at the shore lines. 
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In addition to the dynamic equations (1), the equation of continuity should be included. 
If we neglect the vertical component of velocity, it is 


du dv _ 
dx + dy ” 


(5) 


Equation (5) assumes that there is neither vertical current nor vertical gradient of the 
vertical velocity. Thus, our theory cannot be applied to the coastal and other regions of 
upwelling and sinking caused by local monsoons or other temporary winds. But, if we 
confine ourselves to the gross features of the horizontal circulation in great oceans, induced 
by the superincumbent, quasi-permanent wind system (westerlies or trades), the equation 
of continuity as given by (5) will not cause serious errors in the results. 

There may be some further question concerning the use of (5) for the continuity equa¬ 
tion. But, in treating the oceanographic data for estimating geostrophic currents, we 
always assume 

1 dp 

u = + % t ■ t ; 

Zoo sin <f)p dy 

l dp 

<D = — — 

Zoo sin 4>p dx’ 


provided the frictional terms are neglected. These expressions imply that geostrophic 

currents usually satisfy equation (5) which shows the absence of horizontal divergence. 

^ • • du du » . . 

This means that the equation-1-= 0 may be used without serious errors in treat- 

dx dy 

ing the general circulation of the oceans. 

On eliminating p between the two equations of (1) by cross-differentiation, we have 


. , d 3 u 


d V 
dx 3 




d (du d?A\. 
dz \dy dz/j 


+ 2c op -j- (sin dp) v = 0 
ay 


( 6 ) 


when we take the equation of continuity (5) into account and because — (sin 6) = 0. 

dx 

This equation may also be regarded as expressing the condition that a function p (pressure) 
should exist on a level z as an exact differential with respect to x and y. And the validity 
of equation (6) suggests the possibility of determining the pressure in any level 2 as a 
function of x and y. 


Now suppose the coefficients of mixing Ai and A z are both independent of z, aud put 

V' / \ (25“ — 1) irz 

u = 2s u s (x, y) cos--, (7) 


where 


2 r h 

u s (x, y) = - / u(x, 
‘ Jo 


^ (2s - 1) ^ 

y, f) cos --7^- d£ 


and assume, in accord with Stokes, 


d\ (2s -1) Trf 
5T COS -2h- dt; 


( 8 ) 


( 9 ) 


then we have, by substitution from (2), (3), and (8), 
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/ h 

d 2 u 


COS ^ it - l^ - (2t (X, ,) 


2 h 


h A 2 


4 h 2 


and (9) becomes 

d 2 u _ ^ ^ z t x \x, yj 

ai* " t 3 \h ~T, 

Similarly we have 


/2 r,(x, y) (2s - l) 2 x 2 Ms(X) A cog (2s_- 1) xz 


4h 2 


2 h 


dz 2 


where 


^ (2 r v (x, y) (2s — V/ 7r" / v\ (25 — 1) irz 

f - 5 U - -4P— ’•(*' ( 008 


2& 


V 1 t \ (25 — 1) xz 

W = 2_, V s (x, y) cos--. 


Substitutions of (7), (10), (11), and (12) into (6) and (5) give 


A, 


cfus _ _ (2s - 1) 2 

dy 3 dx 3 


4k 2 


7T 2 A z (du s dvA , ~ cos <f) 

- (57 " to ) + ^ TT ”• 


where jR is the radius of the earth, and 


+ i(r ! -r I) -° 

k \dy dx 


du s . dv 8 = 
dx + dy 

Equation (14) gives a set of functions \p s (x, y) such that 


and (13) becomes 
A /*V , d 4 f.\ (2s - 1) 2 t 2 A z 

Ai kw + w)~~ 


d\p s dips 

U S = -T-; Vs = —— 

dy dx 


4 h 2 


, d 2 \p s \ COS 0 dips 

+ w )~ p ~ ir ~^ 


4.2 /^_ iih = 0 

+ h \dy dx ) ' 


If we introduce two quantities Di and D z such that 

D t — 7 r\l— 

V pQ) 


and 


equation (16) becomes 

-|2 / n 4 , .4 


D z 


x->/— , 

\ pco 


Ei ( d At , §aL] _ (2^ - I ) 2 feY feV. , «V 

2x 2 \ <9x 4 dy 4 / 8 \ h / \ dx 2 dy 2 


cos<£ di/c., 


i? dx 

I 1 ( dry drX _ 

pwA \dy dx / 


( 10 ) 

( 11 ) 

( 12 ) 


(13) 


(14) 


(IS) 


(16) 


(17 


(18) 


( 19 ) 
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The quantity D t may be called “the frictional distance,” whereas D z is the same as 
Ekman’s “depth of frictional influence” except that it does not contain sin <f>. 

The coastal conditions which \f/ s should satisfy are 


Is and = 0 
dn 


( 20 ) 


along the shore lines, where is the derivative of \p s in the direction normal to the 

dn 

shore lines. 

If equation (19) can be solved and we can determine the functions h (x, y), i/^ (x, y), 
fz (x, y), . . . , the sum: 


, / ^ / / \ (2s — 1) 7 tz 

y, z) = 2. I six, y) cos-^- 


( 21 ) 


will give the horizontal streamlines at any level z for the wind stresses r x (x, y) and 
T v y)- The stream function I (x, y, z ) should, of course, satisfy the condition: 


/ d ^ n 
4 ' = Tn = ° 

along the shore lines and the horizontal streamlines of the currents are given by 

p (x, y, z) = constant. 


( 22 ) 

(23) 


Infinitely Deep Ocean 

If i/q is the solution of (19) when h = 1, the solution of (19) will be 

tsix, y) =\ i\(x, y). 

Thus we have the solution (21) of our problem in the form: 

^(x, y, z) = 2^ \ /'Ax, y) ^ cos-^-, 

and we may write down (25) in the form: 

,, \ (x,y) / ttz D Z \ 2 D, 

+ (x ’ y ’ z) = £ “W cos \2^ (2s - 1} TJ —• 


(24) 


(25) 


If we consider the depth h of the ocean increases indefinitely, we have 


( 2 ,— l)fdij (26) 

and 

V (pc, y) —> yp l (x, y; 77 ), (27) 

where fa (x, y, 77 ) is the solution of the equation: 

Ai ( d 4 i/n dVA _ i)_ ( gVi , dVA __ cos 0 dfa J_ / dr* __ dr A _ 

2 cop \dx 4 dy 4 / 8 \dx 2 dy 2 / R dx pw Vdy dx / 

which is derived from (19) by putting (25 — 1) — = 7 ] and h = 1. The right-hand sides 

h 
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of equation (24) will be, when the depth h increases indefinitely, 


fix, y, z ) = 


~ f Hx, y; v) cos (g- v) dr,. 

Z J 0 


(29) 


It may be mentioned that 77 is a parameter increasing from 0 to °°. 


APPLICATION TO THE PACIFIC OCEAN 

In 1950, Munk used the rectangular co-ordinates in discussing the wind-driven oceanic 
circulation in a rectangular ocean. Though he did not take the sphericity of the earth 
into consideration, he could explain the general pattern of the actual circulation of the 
Pacific quite well. Further, his result shows little difference from the one above (Hidaka, 
1951), in which the curvature of the surface of the earth is taken into account and spherical 
co-ordinates are used. 

These results show that the Pacific Ocean can be treated approximately as a rectangular 
ocean, provided we consider that the Coriolis parameter 2 oo sin <j> varies with y. This is 

a quite natural consequence when we consider the relation <£ = where R is the mean 

radius of the earth, and y is counted zero on the equator. 

For these reasons, we can represent the Pacific approximately as a square ocean bounded 

by x ~ 0 , x = a, and y = =±=^ ? in which y = 0 coincides with the equator (Figure 1 ). 

Here a is a mean east-west extent of this ocean and approximately equals 120 ° of longitude 
or 27r/3 in radians. This means that the northern and southern boundaries of this ocean 
2wR 


are given by y = =b— or the parallels of 60° ‘N and S. 
To solve equation (28), assume 

’/' i (%i y, v) = M m {y) N m (x\ vf ), 


(30) 


where m = 1, 2 , 3, . . . and N m (x, rj) are functions of x which are to be determined later, 
while M m (y) are of the forms: 


nr / \ Try miry 

M m (y) = cos-cos — 


for m: odd 


(31) 


ttv . miry 
= COS 2 a Sm ~ 2 a 


for m: even. 

Since these functions and their y-derivatives vanish along y = d=-, (30) and its y-deriva- 


tives will also vanish along y — =j= 5. This makes the function (x, y, z ) satisfy the con¬ 
dition ( 22 ) along the northern and southern boundaries y = ±-. 


Stresses of Winds over the Pacific Ocean 

Dr. Munk kindly furnished me with his unpublished data of the distribution of the 
wind stresses over the Pacific Ocean north of 5° S. From these data, the most probable 
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TABLE I 

Assumed Meridional Variation of the Wind Stress t and Stress Function T 


y/a 

LATITUDE 

T 

dynes/cm 

T 

dyne 1 cm 2 

REMARKS 

yla 

LATITUDE 

T 

dynes / cm 

T 

dyne/ cm 2 

REMARKS 

+ 

.500 

60° N 

.000 



.000 


0.000 

0° 


+ .201 

X 

10 8 

-.272 


+ 

.475 

57° N 

- .037 

X 10 8 

.211 

West 

-0.025 

3° 

S 

.304 

X 

10 8 

-.360 

SE 

+ 

.450 

54° N 

- .141 

X 

10 8 

.396 


-0.050 

6° 

s 

.441 

X 

10 8 

-.463 


+ 

.425 

51° N 

- .301 

X 

10 8 

.528 

Wind 

-0.075 

9° 

s 

.663 

X 

10 8 

-.554 

Trades 

+ 

.400 

48° N 

- .493 

X 

10 s 

.589 


-0.100 

12° 

s 

.811 

X 

10 8 

-.606 


+ 

.375 

45° N 

- .694 

X 

10 8 

.573 

Drift 

-0.125 

15° 

s 

1.018 

X 

10 8 

-.595 


+ 

.350 

42° N 

- .876 

X 

10 8 

.482 


-0.150 

18° 

s 

1.208 

X 

10 8 

-.514 

- . 

+ 

.325 

39° N 

-1.015 

X 

10 8 

.327 


-0.175 

21° 

s 

1.360 

X 

IQ 8 

-.369 


+ 

.300 

36° N 

-1.094 

X 

10 8 

.132 


-0.200 

24° 

s 

1.454 

X 

10 8 

-.174 


+ 

.275 

33° N 

-1.104 

X 

10 8 

-.075 


-0.225 

27° 

s 

1.476 

X 

10 8 

+.045 


+ 

.250 

30° N 

-1.044 

X 

10 8 

-.266 


-0.250 

30° 

s 

1.424 

X 

10 8 

+.259 


+ 

.225 

27° N 

- .926 

X 

10 8 

-.416 

NE 

-0.275 

33° 

s 

1.303 

X 

10 8 

+.438 

West 

+ 

.200 

24° N 

- .767 

X 

10 8 

-.507 


-0.300 

36° 

s 

1.131 

X 

10 8 

+.564 


+ 

.175 

21° N 

- .588 

X 

10 8 

-.534 

Trades 

-0.325 

39° 

s 

.927 

X 

10 8 

+.626 

Wind 

+ 

.150 

18° N 

- .411 

X 

10 8 

-.502 


-0.350 

42° 

s 

.714 

X 

10 8 

+.626 


+ 

.125 

15° N 

- .253 

X 

10 8 

-.431 


-0.375 

45° 

s 

.510 

X 

10 8 

+.574 

Drift 

+ 

.100 

12° N 

- .123 

X 

10 8 

-.343 


-0.400 

48° 

s 

.331 

X 

10 8 

+.484 


+ 

.075 

9° N 

- .024 

X 

10 8 

-.264 


-0.425 

51° 

s 

.187 

X 

10 8 

+.371 


+ 

.050 

6° N 

+ .053 

X 

10 8 

-.219 

Doldrum 

-0.450 

54° 

s 

.083 

X 

10 8 

+.249 


+ 

.025 

3° N 

+ .122 

X 

10 8 

-.219 


-0.475 

57° 

s 

.021 

X 

10 8 

+.124 


+ 

.000 

0° 

+ .201 

X 

10 8 

-.272 


- .500 

60° 

s 

.000 

X 

10 8 

.000 



distribution of the wind stresses was found to be zonal and determined as 

Tx (0) = + 0.045544 M[(4>) ~ 0.262308 Mfo) 

+ 0.022902 ikf'(</>) + 0.069493 M[{4>) (32) 

- 0.036900 Af'(</>) + 0.011560 Mfo) dyne/cm 2 

where $ is latitude so that <f> = y/R and M' m (</>) stand for M m (d>) = R -f- M m (cb). 

. a<t> ay 

If we are to derive (32) from a stress function T (x, y) such that 


dT _ dT 

Tx ~ dy’ Ty ~ dx’ 

we have a stress function which depends on $ or y only, or 


(33) 


T(*) = E T m . M m W 

m =1 

= { 0.29003 Mi(4) - 1.67122 M 2 {4>) 

+ 0.14591 Ms(0) + 0.44276 M^) (34) 

- 0.23510 M & (4>) + 0.07388 M 6 (<t>) } 

X 10 8 dynes/cm, 

provided we take R = 6.3712 X 10 8 cm. 

The values of stress and stress function computed after the formulas (32) and (34) are 
compiled in Table 1 and illustrated in Figure 2. 

From formula (32) we have a much larger wind stress in the west wind belt in the 
South Pacific Ocean than in the North Pacific. It is not because we had wind observations 
available from the former, but is a necessary consequence resulting from determining the 
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formula so as to meet the observational data for the latitudes north of 5° S. As a matter 
of fact, the west wind belt in the South Pacific is said to have stronger wind than that of 
the North Pacific. So this formula will not give wind stresses in the South Pacific too 
inconsistent with observations. 


WEST <- 


-> EAST 


-0.5 


60° N 


45° N 


0.5 


30° S 


45° S 


-0.5 


60° S 


30° N 


15° N 


EQUATOR 


15° S 


WEST WIND DRIFT 


J 


J 

> 


NORTHEAST TRADES 


DOLDRUM 


SOUTHEAST TRADES 


WEST WIND DRIFT 


0 


J 

°- 5 DYNE/CM 2 


Fig. 2. Assumed meridional distribution of wind stress. 

















































192 


PACIFIC SCIENCE, Vol. IX, April, 1955 


Elimination of y-Co-ordinate 

The equation to be solved is from (28) 


pi 

2h 


l\_ (dVi , gV A _ rp (d% dVi\ _ cos 4> dfx , 2 ( dr x dr A = . . 

x 2 \dx 4 dy 4 / 8 \dx 2 dy 2 / R dx poo \ dy dx / 


where 


D l 


- »>/-. 
y poo 


(36) 


Between the equator and 60° N or S, cos $ varies from 1 to 3^. Since this is not a large 
variation, we may take the average value for cos 4> for this range, or 


cos $ 


2 x/3 


cos <j>d(t> 


_ 3_V3 


Thus, the third term in (35) now becomes approximately 


cos <t> d\pi 


3 -\J3 d\f /1 


i? dx 2xi? dx* 
Substituting (30), (32), (34), and (37) in (35), we have 




l(v)‘ 


d 4 N n 


12 7T 2 l(2jr) 4 i? 4 ' dX 4 


M m (<p) + N m (X; 7]) 


1 <f 4 iV„ 


i? 4 * ^. 4 


1 

8 1 ( 25 )' 


^ 2 iV m , Ar A _ , 1 d 2 M, 

(s^) 2 i? 2 ‘ <f\ 2 ' + ’J) ^2 • ^2 




I jit a, \ 1 d Mm 

+ N m (\, v) . R2 . d ^ 2 


d 2 M„ 


3^3 

2 n- 




2_-i?2 * JX 2 

if 




- — — £ T 

pco i? 2 V " d<t> 2 


= 0 


where we have 


X = 


(2irR/3Y 


(37) 


(38) 


(39) 


so that the ocean is supposed to be bounded by the meridians of 0° and 120°, and the 
western and eastern boundaries are given by X = 0 and X = 1 respectively. 

Now it is possible to express ^ and ^ in terms of M m ($), or 


d 2 M n 
dy 2 

d 2 M w 
dy 2 

d A M m 

dy 4 


Substituting (40) in (38) we have 


and 


dy 4 

= £ aTMM) 


= £ jsrAf.-c^)- 


(40) 
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Z Mjfi) 


(2tt 2 ' 




16tr 4 i? 4 d \ 4 


81 


+ Z 


R 2 


v_ 
8 

+ 

^ 2 IT 2 




4tt 2 2? 2 ‘ +y 2 


3 


dN m 

d\ 


• #4 • P™) N i + pu) • £ 2 Z «“ T '} 


= 0. (41) 


In order that this relation should be always valid, the coefficients of M m (<f>) must vanish. 
Thus, we have 

1 DU i N m 


8 _ 

■32 


R 2 ^X 4 


32+ ^ rfX 2 


<7 2 A7„ 


9 V3 r Y' ( 1 i 2 , 1 a i -D?) AT 

“ 1++ v \ ¥ “ m,? + 2+ w Ni 


+ —ZaLT^O (42) 

pco i 

where T* is given by (34). 

It may be anticipated from (39) and (40) that the values of p l m are larger than that of 
a l m for any corresponding set of values of i and m. But the ratio 


B\ = 7 /A 1 

R 2 pco R 2 

is of the order of 10~ 4 for the usual value of Ai (= 10 8 ), so that we may neglect the terms 
of in the coefficients of N{ in (42). Then (42) becomes 


8 1 
■ 52 - 


D\ 


d*N nt 


R 2 d \ 4 


d i N„ 


32tt 2 ' d\ 2 


+ (43) 

4tt 2 d\ Si pco 1 


Determination of the Values of a l m 


The coefficients of in the expansions of 


d 2 M m 

d<t> 2 


the higher terms being neglected. 

For the ocean under consideration, we have 


were evaluated as far as Mi O (0), 


= -itt{Mi(4>) + 8M 3 (4>) - 6M 6 (» + 4M 7 4) - 2 M>(4>)}, 

M 3 '(<3) = —Ai-3M 1 (<t>) + 20M 2 («*>) + 18Af,(0) - 12 M 7 (<f>) + 6M 9 (<3)}, 

Ms(4>) = — t 9 t {+5Ml(4>) - 15M 3 (4>) + 69M 5 (4>) + 20M 7 (<*>) - 10M, (<*>)}, 144) 

M"(cp) = -*{ —7Jlfi(*) + 21M 3 (<3) - 35M 5 (<*>) + 148M 7 (4>) + 14M 9 (<*>)}, 

M 3 (4>) = -*{+9Afifo) - 27 M,(4>) + 45M 5 (<£) - 63M 7 (</>) + 257M 9 (4>)} 

and 

M" (0) = { 7M 2 (0) + 16 M 4 (0) - 12 M 6 (*) + 8I.W-4 M 10 (</>) }, 

M 4 (0) = —Ml —$ M 2 (<3) + 43 M 4 (<£) + 24 Me (4>) — 16 Ms ($) + 8 Mio (4>) }> 

M' t ' (<j>) = — M { +12 M 2 ($) — 24M 4 (4>) + 111 Me (4>) + 24 Ms (4>) — 12 Mio ( <f> ) }, (45) 

M's' (<t>) = — M { —16 M 2 (<*>) + 32 M 4 (<£) - 48 Me (<3) + 211 M 8 (<3) + 16 M 10 (*) }, 

M" (*) = -M { 20 M 2 (</>) - 40 Mi (4>) + 60 M« (0) - 80 Ms (*) + 343 M 10 (0) }. 


Substituting (44) and (45) into equation (43), and evaluating the numerical values of 
their coefficients with R = 6.3712 X IQ 8 cm., pco = 0.000075 sec -1 and 

ill = 3.08 X 10 7 c.g.s., (46) 
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we have for odd series 

D(Ni) + 4747600tj 2 X ttW - 3N 3 + 5 N & - 7 N 7 + 9 N 9 ) 

- 0.548249 X 10 20 = 0, 

D(N S ) + 4747600i? 2 X *(82Vi + 20iV 3 - 152V 6 + 212V 7 - 272V») 

- 3.63163 X 10 20 = 0, 

D(N b ) + 4747600^ X *(-62Vi + 187V 3 + 692V 5 - 35iV 7 + 452V 9 ) (47) 

+ 6.35413 X 10 20 = 0, 

jD(N 7 ) + 4747600 r, 2 X *(+42^ - 122V 3 + 202V 5 + 1482V 7 - 632V 9 ) 

+ 2.19300 X 10 20 = 0, 

D{N t ) + 4747600^ 2 X *(-22Vi + 6N 3 - 107V 6 + 142V 7 + 2572V») 

- 1.096499 X 10 20 = 0, 

where the operator D stands for 

fl-li- 1082323," jl-lWJlloJ;. 48) 

For even series, we have 

D(N t ) + 4747600»j 2 X %(+7N 2 - SN 4 + 12 N 6 - 16N S 

+ 20N 10 ) + 5.45179 X 10 20 = 0, 

D(N 4 ) + 4747600tj 2 X M(+16A7 2 + 437V 4 - 24A7 6 + 322V 8 

- 407Vio) + 3.59829 X 10 20 = 0, 

D(N t ) + 4747600t; 2 X %{-\2N 2 + 242V 4 + 1112V, - 48A7 8 (49) 

+602V 10 ) - 14.76753 X 10 20 = 0, 

D(N g ) + 4747600r; 2 X M(+8 N 2 - 16iV 4 + 24iV 6 + 211A7 8 

- 802Vio) + 7.09512 X 10 20 = 0, 

P(2Vio) + 4747600i) 2 X M(~42V 2 + 8 N 4 - 122V 6 + 16iV 8 

+ 3432V 10 ) - 3.54756 X 10 20 = 0, 

where the operator D is also given by (48). 

To solve the simultaneous differential equations (47) and (48), we employ the method 
of indeterminate multipliers. Let the odd set of equations be of the forms: 

D(N X ) + E C\Ni + Et = 0, 

1=1 

D(N») + E C)Ni + £ 3 = 0, 

1=1 

9 

D(2V,) + E C?2V t - + E, = 0, (50) 

4=1 

D(N 7 ) + E C\Ni + Et = 0, 

4=1 

9 

D(N 9 ) + E C)Ni + £9 = o, 

4=1 

where the summation is made with respect to odd numbers. 
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Multiply each of the five equations in (50) by Z 3 , h , h, and Z 9 , respectively, and add 

together; we have then 


D(l x Ni T I3N3 T kN$ T I7N7 -f" lgN$) 

H - (C1/1 + C1/3 + C1/5 + C1I7 + C1/9) iVi 

+ (Clh + dk + elk + elk + eX) n 3 

+ (clh + elk + elk + elk + elk) n, 

+ (dh + clh + clh + elk + elk) n 7 

+ (dh + clh + elk + elk + elk) n 9 

+ (kEi + kEz + /5-E5 + l 7 E 7 + kET) — 0 


Now let 


elk + elk + elk + elk + elk = he 
elk + elk + elk + elk + elk = he 


elk + elk + elk + elk + elk = he 


and eliminate Zi, Z 3 , . . . , / 9 ; we have then 


Cl-t cl c\ cl cl 

C l /~i3 <_ s~i5 .^t9 

3 C / 3 — £ C 3 C 3 O 3 

/-a a* 5 t /^ 7 

C 5 O 5 O 5 —£ O 5 05 

C 7 1 C 7 3 cl cl -£ C 7 9 

✓~«1 ^-t 9 «. 

0 9 C 9 U 9 U 9 C 9 —£ 


0 . 


(51) 


(52) 


(53) 


This equation has five real roots. Let them be £ 1 , £3 


to them, equations (52) will give five sets of Zi, 

h, . 


£1 

A, 

ii 

k, 

ii 

li, 

£3 

ii 

ii 

i s 

h, 

ii 

li 

£5 

ii 

ii 

ii 

ii 

ii 

£7 

ii 

ii 

ii 

ii 

ii 

£9 

ii 

ii 

ii 

li 

it 


, £5, £7, and £9 and, corresponding 
k, or 


(54) 


These five sets of roots and multipliers l’s will give Y m and F m as 


Vm = IT Nr + IT Ns + IT Ns + IT N 7 + CN 9 ; 

77 _ 7 m Z ? I T m 77 I T? _L T? I 1 m T? 

r m — li JtL\ -j- / 3 rLz -f- h “r h ~r h 
(m = 1 , 3, 5, 7, 9) 


(55) 


and the corresponding five equations 

D(Y m ) + £ m F m + E m = 0 (56) 

(m = 1, 3, 5, 7, 9) 

for Fi, F 3 , . . . , F 9 , They are no longer simultaneous, and are not difficult to solve. The 
same, of course, applies to even sets, too. 
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Practically, the equation corresponding to (53) is 


' 1-f 

+ 8 

-6 

+4 

-2 


-3 

20 —f 

+ 18 

-12 

+ 6 


+5 

-15 

69—r 

+ 20 

-10 

= 0 

-7 

+ 21 

-35 

00 

1 

'OP 

+ 14 


. +9 

-27 

+45 

-63 

257-f 



210830400 - 84156468f + 4625225f 2 - 78771f 3 + 495f 4 - f 6 = 0 (57) 

where 

f = $/(4.38649084 X A) (58) 

and the five roots of (57) are 

fi = 2.9641943, 
f, = 26.7072817, 

f 6 = 74.3718483, (59) 

ti = 146.456904, 
f 9 = 244.499771. 

For even series we have, corresponding to (57), (58), and (59), 


7 — 7 

+ 16 

-12 

+8 

-4 


-8 

43- 7 

+ 24 

-16 

+8 


+ 12 

-24 

111-7 

+ 24 

-12 

= 0 

-16 

+ 32 

-48 

211 —y 

+ 16 


+ 20 

-40 

+ 60 

-80 

343- 7 



5217079023 - 593115237 7 + 16671798+ - 172458+ + 715+ - 7 5 = 0 (60) 

where 

7 = £/(4.38649084 X %) (61) 

and 

72 = 12.8567852, 

74 = 51.494636, 

7 6 = 116.1536052, (62) 

7s = 207.4231320, 

7 io = 327.071842. 


Differential Equations for Y's and Their Solutions 

i 

From the numerical computations described above, the numerical coefficients of the 
differential equations for F m (A; rj) were determined as follows: 

D(Fi) + 11514116 j) 2 Fi - 0.4037686 X 10 20 = 0, 

D(Y 3 ) + 103741 765ij 2 F 3 - 6.3254121 X 10 20 = 0, 

D(Y b ) + 288890008 r, 2 Y b - 204.1449265 X 10 20 = 0, (63) 

£>(F 7 ) + 568897198i7 2 F 7 + 311.4127463 X 10 20 = 0, 

£>(F 9 ) + 949734910»; 2 F 9 + 257.7413034 X 10 20 = 0 
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D(Y t ) + 51883042r/ 2 V 2 + 7,4929828 X 10 20 = 0, 

D(Y 4 ) + 207804544+ F, - 5.3393918 X 10 20 = 0, 

D(Y 6 ) + 468733228^ F 6 + 258.843806 X 10 20 = 0, (64) 

D(Ys) + 837047752+F 8 + 158.8661390 X 10 20 = 0, 

D( Fio) + 1319885336+Fio + 290.043938 X 10 20 = 0 


where 

j2 7 

D = —7 - 1082323+ -r-r - 14997110-,- 
dX 4 ' d\ 2 d\ 

and 


Fi, Y s , , F 9 ; F 2 , Y 4 , ... , Fio stand for the following expressions in terms of 7V m : 
F x = N 4 + 0.21670317^3 - 0.0302211TV 5 + 0.0100338TV 7 - 0.00455397V 9 , 

F 8 = 7V X + 4.63327077V 3 + 1.67781057V S - 0.2720794TV 7 + 0.10185167V 9 , 

Y*> = N\ — 7.3461297^3 - 31.5044500TV 5 - 13.1812666TV 7 + 2.08037397V 9 , (65) 

Yt = Ni — 4.28834687V 3 + 20.40048077V 5 + 99.04981447V 7 + 47.01634757V 9 , 

F 9 = 7V X - 3.65835477V 3 + 9.78495227V 6 - 40.09696637V 7 - 245.93209357V 9 
d 

F 2 = 7V 2 + 0.31722317V 4 — 0.0501698TV 6 + 0.0180624TV 8 - 0.0086497TVi 0 , 

F 4 = 7V 2 + 4.12626137V, + 1.58741777V 6 - 0.25978617V 8 + 0.09956097Vio, 

F 6 = 7V 2 - 5.4990870TV, - 23.8405548TV 6 - 10.31299497V S + 1.61092537Vi 0 , (66) 

F 8 = 7V 2 - 3.08742347V, + 13.82854477V 6 + 68.62958777V S + 33.31806477V 10 , 

Fio = 7V 2 - 2.57126547V, + 6.42713777V 6 - 25.69390467V S - 160.97224447V 10 . 


If the 10 functions Fi, F 2 , . . . , Fi 0 can be determined by solving the differential equa¬ 
tions (63) and (64), it will be possible to compute 7Vi(X; ij), 7V 2 (X; »;),..., Nio(X; ?j) from 
the following expressions which are the reversions of the expressions (65) and (66). 


and 


TVi - +1.0650221 Fi - 0.05878696F 3 - 0.00491292F 5 

- 0.001038174F 7 - 0.000284101 F 9 , 
7V 3 = -0.2793058 Fi + 0.2586992 F 3 + 0.01644431 F 5 

+ 0.003283144F? + 0.000879072 F*. 
7V 5 = +0.1254434F! - 0.07711532 F 3 - 0.04050495 F 6 

- 0.006251173 F 7 - 0.001571970 F 9 , 
TV 7 = -0.05969848 Fj + 0.03367969 F 3 + 0.01083299 F 6 

+ 0.01265900 F 7 + 0.002526788 F 9 , 
7V 9 = +0.02320969 Fi - 0.01264167 F 3 - 0.003642393 F 5 

- 0.002365708 F 7 - 0.004554906 F 9 


(67) 


7V 2 = +1.1172140F 2 - 0.1028957F, - 0.0109283F 6 

- 0.0026156 F 8 - 0.0007745 Fio, 

TV, = -0.3385977 F 2 + 0.3066091 F, + 0.0247985 F 6 

+ 0.0055792 F 8 + 0.0016107 F w , 

TVo = +0.1613919 F 2 - 0.0943220 F, - 0.0548464 F 6 

- 0.0096171 F 8 - 0.0026064 Fio, (68) 

7V 8 = -0.0792995 F 2 + 0.0420747 F, + 0.0147581 F 6 

+ 0.0184664 F 8 + 0.0040001 Fio, 

TVio = +0.0314504 F 2 - 0.0160186 F, - 0.0050095 F 6 

- 0.0034369 F 8 - 0.0069853 Fi 0 . 
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Substituting the functions Ni(x; 77 ), N 2 (x; 77 ), , Ni 0 (x; 77 ) thus obtained in (30), we 

shall obtain the expression for yp\(x, y; rj). Further substitution in (29) will give the solution 
of the problem as 

i{x, y, z ) = M m {y) . j N m {x\ rf) cos (ff- dr). (69) 

The method of solving equations (63) and (64) and the evaluation of the integral: 


r 


N m (x; rf) cos 




drj 


will be discussed in the following section. 


Integration of the Differential Equations 


The next step will be to solve the 10 differential equations (63) and (64). Let any one 
of these equations be 

d*Y 
d \ 4 


1082323j? 2 

dh z 


dV 

14997110 -j- 
a\ 


+ brj Y + c = 0 


(70) 


where b and c are constants assigned to each of these 10 equations. Since a particular 
solution of this equation is 


Y p = 



(71) 


the general integral of (70) will be of the form: 

F(X) = Ac aX + Be^ + Cc yX + Dc SX - ( 72 ) 

where a , jS, 7 , and 5 are the four roots of the algebraic equation: 

<r 4 - 108232377V - 149971 lOer + brj 2 = 0 (73) 

and A, B, C, D are constants to be determined according to the conditions: 

F( 0 ) = 7(1) = 0 ; Y[( 0 ) = Y[( 1) = 0 . (74) 


The equation (73) has four roots a , 7 , and 5 for any given value of 77 , and A, B, C, 

and D all depend upon a , (3, 7 , and 3. 

The parameter 77 varies from 0 to 00 , and the values of a , 7 , 8 all depend upon 77 . 

When 77 is very large, these four roots are approximately 

a = +1040.307577, 

P = -1040.307577, (75) 

7 = +(6/1082323)*, 

8 = -(6/1082323)1 

As 77 decreases, a , ft, 7 , and 8 also change gradually. For 77 less than a certain value be¬ 
tween 77 = 0.4 and 77 = 0.3, @ and 8 become complex conjugate. As 77 approaches 0, a , 
ft, and 8 approach finite values, while 7 decreases indefinitely as «: +. Thus we have, 
when 77 —> 0 , 

a = 2 p, 

P = —p + iq , 

7 = brj 2 /Sp z , 

8 = —p — iq 


( 76 ) 
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where 

p = (14997110)*, 



(77) 


For the intermediate value of % 7 , these four roots vary continuously except 0 and 8 which 
change from complex conjugates to real as rj increases from rj = 0.3 to r} = 0.4. 

Of course, there are 10 series of such four roots of 10 equations (70), each varying with 
the parameter 17 . 

The constants A , B f C, and D can be, of course, expressed in terms of a, 0 , 7 , and 8, 
Thus the solution F becomes, when fi and 8 are complex conjugates of the form 


¥ = jl + —e“ (X_1> -— 

l a — 7 a — y 


P = +P + qi; 5 = p — qi, 

oc 7(X_l) / a -7 , \ — pX 

+ ^ ■ _ ■ ■ e — l J e cos qk 


(78) 


and, when a, j 8 , y, S are all real, 
F = 


a — 7 


+ (—■ p ~y - *) 4 • r» (79) 

}• 


! + - 2 — e“ (X-l) + ^ ^ e~ 7 } . / X 


.0 — 8 a — y * fi — 8 


a t(x_i) . 

- e + 

OL — 7 


0 0-7 


/5 — d a — y * fi — 8 


px 


. e 


where F p is the particular solution given by (71), of the equation (70). 

When 7] increases, 7 also increases. If we can neglect e ~ T , the expression (80) will be 
further simplified, and we have 


Y= il + 


7 a(X_i) 5 

-e + 


« — 7 fi — 5 

Since Fp is given by (71) as 


/3X 

6 — 


« 7(X_i) 

- g 

a — 7 


0 sx 
£ 


0 “ 5 


} • ¥,. 


(81) 


F» = — 


(c/ 6 ) 


the solutions of the equations (63) and (64) will tend to zero as rj increases indefinitely. 

The values of the roots a , fi , 7 , and 8 of each of the 10 equations given by (63) and 
(64) were computed numerically and are given in Tables 5-14. 


Computation of the Currents 

In order to calculate the distribution of the currents at various levels, we had first to 
compute Fi(A; 77 ), F 2 (A; 77 ), ... , F 10 (A; rj) according to one of the expressions (79), (80), 
and (81) for 

A - 0.0000, 0.0025, 0.0050, „ . . , 0.0500, 0.0550, 0.0600, . . . , 0.1000, 0.1100, 

0,1200, 0.1300, . . . , 0.2000 

and for 

77 - 0.0, 0.1, 0.2, . . . , 1.0, 2.0, 3.0, . . . , 10.0. 
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For larger values of r], we have 

a = + 1040.3075*?, 

/3 = — 1040.3075*7, (82) 

y = + (&/ 10S2323) 

8 = - (6/1082323) 1 

very accurately, while Y p (r 7 ) are all very small, so that F 2 (A; 77 ) will be approximately 
given by 

Y = jl — e T(X_l) - e _7X } . (83) 

where 7 , and c are independent of 77 . It will be more convenient to leave (83) as it stands 
rather than to compute their values against rj. 

These values of the 10 functions F m (A; 77 ) may be, then, converted into the functions 
N m (\; r]) by virtue of formulas (67) and ( 68 ). 


Substitutions of the functions iVi(A; 77 ), iV 2 (X; 77 ), iV 3 (X; 77 ), ... , iVi 0 (X; 77 ) into equation 
(69) give the complete solution as 


y, z) = 

To evaluate the integral 



(84) 


we have computed the functions 2V m (A; 77 ) for 

77 = 0.0, 0.1, 0.2, . . . , 1.0, 2.0, 3.0, . . . , 10.0 
and the process has to be carried out numerically between 77 = 0 and 7 ? = 10 . 0 . For 
larger values of 77 , we may use the approximate formula (83). 


Let the values of a function F( 77 ) for 77 = 0 , h, and 2 h be F 0 , Fi , and F 2 , respectively. 
Then the interpolation formula in this interval of 77 will be given by 


F( 7?) = F 0 + 
Then we have 

^2 h 


— 3F 0 T 4 Fi — F 2 ( 77 \ ^ Fp — 2Fi T F 2 [ 77 


/ 


or 


F( 77 ) cos 


dr) + 


h/ 


— 3Fq + 4Fi — F 2 


(85) 


fe v ) dv = Fo f cos fe v ) 

•Sq 

r 2h ( 7TS ^ J , Fo - 2F, + F 2 r 2h 2 _ ( *Z \ A 

/ y ' COS \2D 1 V dr ' + -~~W -/ " COS \2Dj r,dri 

•So •So 


f F(rj) cos vj dr) = (ai — 3bi + £ 1 ) • F 0 + (4&i — 2ci) . Fi + (&i + £ 1 .) • F 2 
•So z 
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where 


ai 


-J cos (w, ’) 


= —-— sin (—■ . 2k) 

/jr z_\ \2D Z J 

\2D Z ) 

bi ~2hf vcos (m v ) dv 

•SQ 

1 (it ■ 2 *) ,in (it • 2h ) + cos (it • 2h ) - 1 

(' m ;) (25; • 2 *) 

■4 s*2 h 

Cl ~ 2h?f ' COS (W z V ) drl 

1 2 (it ■ 2k ) ™ (it ■ 2h ) + \(w ,- 2h )~ 2 \ sin (it ■ 2h ) 


( 86 ) 


(it) 


(it ■ 2 ‘) 


For larger values of rj , we may express in the form 


F(v) = ft + ~ F ~‘ (, - w) + r+ ' ~ T r ~‘ (, - „)* 

where F_ 1 , F 0 , and F +i are the values of F(^) for rj 0 — h } rj 0 , and ry 0 + A, respectivelv. 


F + i — 2Fq + F 


We have then 


/ '/OT" 

-F(t?) cos (~ 2 ]j~ ) ~ (“~Z ?2 + c 2 ) F_i + (a 2 — 2c 2 ) F 0 + (& 2 + c 2 ) F+i 

.-h 


1871 


where 


a 2 = 11 — I ( 

> irz 

<2 17/ 

*)1 

cos ^ 

' 7TZ \ 

JD Z V ‘ 

2 h, 


s - - i* (it 

■*)• 

” ¥ 1 o 1 

/ 7TZ 

K2D, 

. hy\ sin 

(w. *>) h ’ 

(88) 

C 2 = | f “ TO 

/ TZ 

\2Dl 

• 

| cos 


. h. 



The integration was carried out taking h = 0.1 in the interval 0 S y ^ 1.0 and h — 1 
for the interval 1.0 g rj S. 10.0. For larger values of rj, integration was carried out by 
using formula (82) which is only inversely proportional to rj 2 . In this case, we have only 
to evaluate the integral: 




( 89 ) 
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TABLE 2 

Computed Velocity* of the Kuroshio at Different Depths z = 0, V 2 D z , D z , 3 / 2 D z , 2D Zi and 3D 2 
Along the 33° N Parallel, Assuming D z = 75 meters 


Distance from west boundary (km.) . . . . 

. . . 0 

25 

50 

75 

100 

125 

150 

175 

200 225 250 

275 

Surface velocity (cm/sec). 

. . . 0 

163 

217 

196 

146 

94 

52 

25 

12 7 8 

11 

Velocity at z = Vi D z (cm/sec). 

. . . 0 

150 

206 

186 

135 

84 

43 

17 

4 0 1 

4 

Velocity at z = D z (cm/sec). 

. . . 0 

138 

188 

167 

120 

72 

32 

7 

-5 -8 -6 

-3 

Velocity at z = 3 / 2 D z (cm/sec). 


123 

168 

149 

104 

47 

21 

9 

-12 -14 -12 

-8 

Velocity at 2 : = 2 D z (cm/sec). 

. . . 0 

107 

147 

129 

87 

44 

11 

-9 

-22 -19 -15 

-11 

Velocity at z = 3 D z (cm/sec). 

. . . 0 

79 

110 

94 

59 

24 

-3 

-25 

-25 -17 -18 

-13 


* To get the velocity for values of D z other than 75 m. multiply these velocities by 75 /D z , where D z is ex¬ 
pressed in meters. 


The computations were rather tedious and took three computers more than 6 months to 
complete for the surface, D Z1 %D Z , 2 D z , and 3D Z . The values of a, 7 , and d , which 

are the roots of the quadratic equations (73) with rj as a parameter and computed for the 
necessary values of rj, are compiled in Tables 5-14. 

The values of the stream-function were computed for the westernmost one-fifth part 
of the entire expanse of the ocean, and the streamlines were drawn for the layers z = 0 
(surface), x /iD z , D z , \D Z , 2D Z} and 3D Z . The computations were not carried out for the 
deeper levels and for the part to the east of this area, partly because we did not have 
enough time to compute, and partly because the central and eastern part are not as 
interesting. We have only a very slow zonal flow in the central part and very diffuse 
meridional flow close to the eastern coast of the ocean. Indeed, the California and Peru 
Currents are considered to be produced by local winds as proved by Munk (1950). 

The circulation patterns in the area close to the western coast were obtained from these 
computations and illustrated in Figures 3a-f. The discussions for them will be given in 
the following paragraphs. 

Table 2 gives the velocities of the western current in the subtropic gyre corresponding 
to the Kuroshio, or the Japan Current, at the depths .z = 0 (surface), y%D z , D Zl \D Z , 
2 D ZJ and 3D Z along the 33° N parallel which is the swiftest part of this mighty current. 
These velocities were computed by the formula: 


v = 



^(Ao T AX) — yp (Ao — AX) 
linear distance of 2AX 


and taking D z — 75 meters. Because the velocity is inversely proportional to the quantity 
D Z1 we can compute it for any other value of D z . For this we have only to multiply these 
figures by 75 /D z where D z is expressed in meters. The maximum surface velocity of 217 
cm/sec agrees with actual observations very closely. 

Table 3 gives the distribution of E-W components along a meridian 24° of longitude 
east of the western boundary, or one fifth of the entire east-west expanse of the Pacific 
Ocean off the western coast. At this distance from the western coast, the coastal effect 
nearly vanishes and the pattern of the circulation consists of approximately zonal flows. 
In this table the value of D z was again assumed to be 75 meters. Discussions concerning 
these results will be made in the following paragraphs. 

Surface Circulation 

The numerical result for the horizontal circulation has been worked out for several 
levels specified by the ratio z/D z . We show here those of the surface (z = 0), z/D z = 0.5, 
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z/D z — 1, z/D z — 1.5, z/D z — 2, and z/D z = 3. The most important of them is, of 
course, the surface circulation, and Figure 3 a shows its pattern. The gross features of the 
current distribution on the surface thus do not seem to differ much from those given by 
Munk (1950) and by the author (1951) for the distribution of mass transport streamlines. 


TABLE 3 

Zonal Distribution of E-W Components of Ocean Currents Expressed in cm/sec 
Assuming D z = 75 Meters (+ Eastward, — Westward) 


LATITUDE 

z = 0 
(surface) 

Z = !/2 D z 

z = D z 

z = 3 A d 2 

z = 2 D z 

z = 3 D z 

REMARKS 

60° N 

0 

0 

0 

0 

0 

0 

Arctic 

57° N 

0 

- 1.2 

- 2.1 

- 2.7 

- 2.9 

- 2.8 

Current 

54° N 

+ 2.5 

+ 0.2 

- 1.7 

- 2.8 

- 3.5 

- 3.6 


51° N 

+ 8.2 

+ 5.0 

+ 2.3 

+ 0.5 

- 0.7 

- 1.7 


48° N 

+ 15.3 

+ 11.6 

+ 3.6 

+ 6.1 

+ 4.2 

+ 2.3 

West 

45° N 

+20.4 

+ 16.7 

+ 13.3 

+ 10.7 

+ 8.7 

+ 5.9 

Wind 

42° N 

+20.6 

+ 17.4 

+ 14.5 

+ 12.2 

+ 10.3 

+ 7.3 

Drift 

39° N 

+ 14.9 

+12.8 

+10.8 

+ 9.2 

+ 7.8 

+ 5.7 


36° N 

+ 5.0 

+ 4.3 

+ 3.6 

+ 3.0 

+ 2.6 

+ 2.0 


33° N 

- 5.8 

- 5.0 

- 4.4 

- 3.8 

- 3.3 

- 2.3 


30° N 

-14.1 

-12.1 

-10.3 

- 8.8 

- 7.4 

- 5.3 

North 

27° N 

-18.3 

-15.4 

-12.8 

-10.7 

- 8.8 

- 6.0 

Equatorial 

24° N 

-18.4 

-15.1 

-12.2 

- 9.8 

- 7.9 

- 5.1 

Current 

21° N 

-15.8 

-12.5 

- 9.8 

- 7.5 

- 5.8 

- 3.5 


18° N 

-11.8 

- 8.9 

- 6.5 

- 4.6 

- 3.3 

- 1.7 


15° N 

- 7.3 

- 5.0 

- 3.2 

- 1.9 

- 1.1 

- 0.4 


12° N 

- 2.7 

- 1.2 

0 

+ 0.7 

+ 1.0 

+ 0.8 


9° N 

+ 2.4 

+ 3.3 

+ 3.9 

+ 4.0 

+ 3.8 

+ 2.8 

Equatorial 

6° N 

+ 7.5 

+ 7.0 

+ 9.3 

+ 7.1 

+ 6.5 

+ 4.8 

Counter 

3° N 

+ 7.6 

+ 8.2 

+ 8.4 

+ 8.2 

+ 7.5 

+ 5.7 

Current 

0° 

+ 5.3 

+ 6.3 

+ 7.0 

+ 7.0 

+ 6.7 

+ 5.2 


3° S 

- 0.7 

+ 1.6 

+ 2.5 

+ 3.2 

+ 3.5 

+ 3.1 


6° S 

- 6.9 

- 6.4 

- 4.1 

- 2.5 

- 1.5 

- 0.5 


9° S 

-16.4 

-12.9 

- 9.8 

- 7.5 

- 5.8 

- 3.5 

South 

12° S 

-22.9 

-16.8 

-13.3 

- 4.9 

- 8.3 

- 5.2 

Equatorial 

15° S 

-21.5 

-17.6 

-14.1 

-11.4 

- 9.1 

- 5.9 

Current 

18° S 

-18.3 

-14.9 

-11.9 

- 9.5 

- 7.6 

- 4.8 


21° S 

-12.6 

-10.2 

- 8.0 

- 6.3 

- 5.0 

- 3.0 


24° S 

- 5.9 

- 4.8 

- 3.7 

+ 3.0 

- 2.3 

- 1.4 


27° S 

+ 1.2 

+ 0.8 

+ 0.6 

+ 0.4 

+ 0.3 

0 


30° S 

+ 8.2 

+ 6.5 

+ 5.0 

+ 3.9 

+ 3.0 

+ 1.7 


33° S 

+ 14.9 

+ 12.0 

+ 9.5 

+ 7.5 

+ 5.8 

+ 3.5 

West 

36° S 

+20.0 

+ 16.5 

+ 13.2 

+ 10.6 

+ 8.4 

+ 5.3 

Wind 

39° S 

+22.4 

+ 18.4 

+ 14.8 

+ 12.0 

+ 9.7 

+ 6.3 

Drift 

42° S 

+20.0 

+16.4 

+ 13.1 

+ 12.5 

+ 8.4 

+ 5.6 


45° S 

+13.7 

+ 10.6 

+ 7.9 

+ 6.0 

+ 4.7 

+ 3.0 


48° S 

+ 5.0 

+ 2.8 

+ 1.0 

0 

- 0.5 

- 0.6 


51° S 

- 2.6 

- 3.4 

- 4.8 

- 5.1 

- 4.9 

- 3.8 

Antarctic 

54° S 

- 6.2 

- 6.8 

- 7.2 

- 8.0 

- 6.4 

- 4.8 

Current 

57° S 

- 4.8 

- 5.0 

- 5.1 

- 4.9 

- 4.4 

- 3.3 


60° S 

0 

0 

0 

0 

0 

0 
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Because of the considerable labor contained in the calculation, the computation is con¬ 
fined only to the western part of the ocean bounded by two meridians, X = 0 and X = 0.2, 
that is, 24° of longitude apart. Choice of the western part of the ocean for the computation 
was made because the circulation patterns in that section are more complicated and hence 
more interesting. In the central part of the ocean we will have indeed a very slow motion 
approximately in east-west direction, while very diffuse meridional motion will exist close 
to the eastern coast. 

We have a number of gyres in the surface circulation corresponding to those obtained 
in Munk’s (1950) and the author’s (1951) results with respect to mass transport. We have 
a broad gyre with strong western current flowing north in the latitudes between 20° and 
40° N and corresponding to the Kuroshio, or Japan Current. We also notice one boundary 
vortex, but the secondary boundary vortex is not distinct. We have a subtropic gyre with 
the western current flowing south. This corresponds to the Mindanao Current. Of course, 
we have a faint subarctic gyre corresponding to the Oyashio, or the Kurile Current. 

On the surface of the Southern Pacific Ocean, we have western currents flowing north 
a little south of the equator and in the subantarctic latitudes. Between these two we have 
a strong current corresponding to the East Australian Current, though actually this 
current never develops so strongly because of many passages connecting the Southern 
Pacific to the Indian Ocean through the numerous islands and archipelagoes in the 
Australian-Asiatic Mediterranean. Had we not these passages together with the Southern 
Antarctic Circumpolar Ocean, we could have a much stronger western current in the 
South Pacific Ocean than actually observed. 

It looks also rather strange that we do not have any strong westward flow in the 
latitudes between 5° N and 2° S although actually the northern margin of the South 
Equatorial Current is in this zone. This is because the Equatorial Counter Current appears 
in our theoretical result much broader and much more diffuse than actually observed. 
This is also the same in Munk’s and Hidaka’s results. The theory of the Equatorial Counter 
Current has been attacked and explained by several authors (Montgomery and Palmen, 
1940; Neumann, 1947) in some other ways than ours. 

Evaluation of the Coefficient of Vertical Mixing 

The streamlines in Figure 3 a are drawn for an interval A\(/ = 250 X 10 10 /D Z cm 2 /sec 
of the stream function. The velocity can be determined as the ratio A\p/Ax, where Ax is 
the actual distance between two consecutive streamlines. Because these diagrams are not 
drawn in a common scale for the north-south and east-west directions, it would be 
rather laborious to compute the magnitude of current velocity for all parts of the Pacific. 
Still it will be easy to determine when the streamlines run in exactly north-south or 
east-west directions. 

The values of the stream function at several points along the 33° N parallel are com¬ 
puted as compiled in Table 4. Assuming the Pacific Ocean is 10,000 kilometers across in 
its east-west direction, we obtained the velocity of the Kuroshio at its swiftest zone, which 
is located approximately 55 km. off the coast, to be 329 cm/sec, 219 cm/sec, 165 cm/sec, 
and 110 cm/sec according as we assume D z = 50 m., 75 m., 100 m., and 150 m., respectively. 

Actual velocity of the Kuroshio has been estimated at approximately 3 to 5 knots, or 
about 150 to 250 cm/sec in its swiftest zone. From Table 4 we recognize that the computed 
velocity of the Kuroshio, assuming for D z a value between 50 m. and 150 m., agrees with 
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Fig. 3. Computed circulation in the western part of the Pacific Ocean at six levels, it. Surface; b, z = Vi D z \ c, 2= 1.0 D-; d, z= 1.5 D,\ e, z=2 D : \ f, 2 = 3.0 D-. 
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TABLE 4 

Computation of Theoretical Surface Velocity of the Kuroshio across the 33° N Parallel 


X = 

Distance from 

west boundary (km.) = 

0.0025.0050.0075.0100.0125 .0150 .0175 .0200 .0225 .0250 .0275 

0 25 50 75 100 125 150 175 200 225 250 275 



* = 

0 

612 

815 

735 

546 

352 

195 

95 

44 

28 

31 

41 

X10 10 cm 2 /sec 















dT 


A\p/Ax = 

0 122.4 163.0 157.0 109.2 

70.4 

39.0 

19-0 

8.8 

5.6 

6.2 

8.2 

X10 4 cm/sec 















dT 

Computed , 

(For D z = 50 m. 

0 

245 

326 

314 

218 

141 

78 

38 

10 

11 

12 

16 


Velocity , 

jFor D z = 75 m. 

0 

163 

217 

207 

146 

94 

52 

25 

12 

7 

8 

11 


(cm/sec) 

| For D z = 100 m. 

0 

122 

163 

157 

109 

70 

39 

19 

9 

6 

6 

8 



' For D z = 150 m. 

0 

82 

109 

105 

73 

37 

26 

13 

6 

4 

4 

5 



the observed values fairly well. The previously determined values for D z fall mostly in 
this range also. This enables us to compute the values of vertical coefficient of mixing from 
formula (18). The above values of D z correspond to the values 188, 422, 750, and 1688 
g/cm/sec of A Zf respectively. These are, of course, values consistent with the results de¬ 
rived from many other different sources. (Sverdrup, et al., 1942.) 


Subsurface Circulation 

Figure 3b, g, d, e, and f show the horizontal distribution of streamlines in the level } [ /^D z , 
D z , 1 }/ 2 D z1 2 D z , and 3D Z below the sea surface, respectively. All give patterns similar 
to the sea surface circulation shown in Figure 3a. We have western currents and a boundary 
vortex attached to each gyre. The only difference noticed is a general subsidence of the 
motion as we go down into deeper layers. Still, we see that the intensity of motion is only 
reduced to as low as half that of the sea surface even in the layer 3D Z . Figure 4 shows the 
comparison of the current velocity profiles along the 33° N parallel at several levels to 
that on the surface of the sea assuming D z = 75 m. The maximum speeds are seen at 
about 55 kilometers off the western boundary. Although the Japanese Islands are not 
disposed parallel to a meridian, the above result agrees with the observed profiles of this 
mighty current quite satisfactorily. Another result of particular interest is that, at a 
distance greater than about 150 km., there is a flow to the south with much greater velocity 
than in upper layers. This counter current reaches a maximum speed of 20 to 30 cm/sec 
at about 200 km. off the western coast, despite the practically motionless upper layers. 
Figure 5 gives the comparison of the zonal distribution of E-W components along a merid¬ 
ian 24 degrees of longitude to the east off the western boundary. In this longitude it is 
expected that the influence of the western boundary nearly vanishes and the actual flow 
pattern of the Pacific circulation is disposed mostly as a zonal current system. The vel¬ 
ocities of the current in these diagrams were computed assuming D z — 75 meters. For 
computing the velocities when the value of D z is different, we have only to multiply these 
figures by 75 /D z , where D z is expressed in meters. If we assume, however, that the value 
D z — 75 m. is consistent, we have for the maximum surface velocities of the North Pacific 
Current, North Equatorial Current, Equatorial Counter Current, South Equatorial Cur¬ 
rent, and Antarctic Circumpolar Current 22, 19, 8, 23, and 23 cm/sec, respectively. They 
are reduced to 18, 16, 8, 18, and 18 cm/sec, respectively, at a level }/ 2 D z and to 11, 9, 8, 
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Fig. 5. Zonal distribution of E-W components of ocean currents at several depths expressed in cm/sec, when 
D z — 75 meters. (To obtain the velocity for other values of D z , multiply these velocities by 15/D z , where D z is 
expressed in meters.) 


9, and 10, respectively, at 2 D z , The Equatorial Counter Current remains nearly unaltered 
in its speed in all depths compared above. 

Vertical Variation of the Currents 

The most important objective of the present research is to get a certain idea about the 
vertical structure of the wind-driven circulation in the Pacific Ocean. This will be, of 
course, impossible to obtain from former theories which have been mostly propounded 
with respect to the mass transport. 

As the problem is three-dimensional, the numerical computation is rather tedious. For 
this reason the author has not yet been able to finish the computation below the level 
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z = 3D Z . Still, we should be able to expect some important conclusions from what has 
been completed thus far. 

First of all it is very interesting that the wind-driven currents exist in a layer much 
deeper than that expected from the classical theory of Ekman (1905). According to Ek- 
man’s theory, a wind-driven current is confined to the surface layer about D z thick, and 
we can expect practically no drift current at a deeper level than 150 m. except very close 
to the equator. From our computation, it can be shown that the current velocity does not 
drop as low as half the surface value even at a level 3D Z . If we take D z = 75 meters, this 
depth is 225 meters. 

This conclusion will help us understand the fact that wind-driven currents can pene¬ 
trate into a layer several hundred meters deep, several times as deep as that expected 


TABLE 5 

Four Roots of the Equation: <r 4 — 1082323??V 2 - l4997li0«r + 1151411 6+ = 0 (m = 1) 
For Different Values of r\. Refer to Equation (73) 


V 

a 

0 

y 


8 

0.1 

261.21657 < 

1 - 130.612124 

\ - 200.881495 i 

0.0076775 

S -130.612124 
/ +200.8814951 

0.2 

304.26551 \ 

> - 152.148109 

\ - 161.6954191 

0.0307072 

i -152.148109 
( +161.6954191 

0.3 

371.21622 \ 

j — 185.642643 

j- 77.1685151 

0.0690668 

i; 

185.642643 

77.1685151 

0.4 

454.08159 

- 363.12143 

0.122667 

— 

91.082828 

0.5 

545.93222 

- 489.85023 

0.191279 

— 

56.273270 

0.6 

642.62244 

- 603.98223 

0.274436 

— 

38.914645 

0.7 

741.98396 

- 713.66160 

0.371324 

— 

28.693687 

0.8 

842.89309 

- 821.22710 

0.480691 

— 

22.146678 

0.9 

944.74621 

- 927.63340 

0.600783 

— 

17.713594 

1.0 

1047.20274 

- 1033.34356 

0.72937 

— 

14.58855 

1.1 

1150.06112 

- 1138.60821 

0.86382 

— 

12.31673 

1.2 

1253.19657 

- 1243.5733 6 

1.00136 

— 

10.62457 

R3 

1356.52895 

- 1348.32950 

1.13922 

- 

9.33867 

1.4 

1460.00495 

- 1452.93512 

1.27490 

— 

8.34473 

1.5 

1563.5883 

- 1557.4297 

1.40631 

— 

7.56490 

1.6 

1667.2526 

- 1661.8399 

1.53190 

— 

6.94465 

1.7 

1770.9803 

- 1766.1856 

1.65059 

— 

6.44526 

1.8 

1874.7576 

- 1870.4809 

1.76177 

— 

6.03849 

1.9 

1978.5740 

- 1974.7356 

1.86521 

— 

5.70361 

2.0 

2082.4224 

- 2078.9583 

1.96096 

— 

5.42509 

2.5 

2601.9746 

- 2599.7576 

2.33637 

— 

4.55339 

3.0 

3121.8104 

- 3120.2708 

2.58146 

— 

4.12107 

3.5 

3641.7804 

- 3640.6492 

2.74475 

— 

3.87598 

4.0 

4161.8219 

- 4160.9559 

2.85725 

— 

3.72328 

4.5 

4681.9050 

- 4681.2207 

2.93742 

— 

3.62168 

5.0 

5202.0138 

- 5201.4595 

2.99627 

— 

3.55053 

5.5 

5722.1396 

- 5721.6815 

3.04065 

— 

3.49871 

6.0 

6242.2768 

- 6241.8919 

3.07487 

— 

3.45977 

6.5 

6762.4223 

- 6762.0943 

3.10178 

— 

3.42975 

7.0 

7282.5734 

- 7282.2906 

3.12332 

— 

3.40610 

7.5 

7802.7290 

- 7802.4827 

3.14080 

— 

3.38714 

8.0 

8322.8879 

- 8322.6714 

3.15518 

— 

3.37170 

8.5 

8843.0494 

- 8842.8577 

3.16716 

— 

3.35890 

9.0 

9363.2128 

- 9363.0417 

3.17724 

— 

3.34830 

9.5 

9883.3778 

— 9883.2242 

3.18579 

— 

3.33937 

10.0 

10403.544 

-10403.405 

3.19310 

— 

3.33166 
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from Ekman’s theory of the limit of the wind-driven currents. This implies that the mo¬ 
tion of water in most parts of the oceanic troposphere could be produced by the stresses 
of the permanent wind system prevailing over the oceans. In other words, the winds are 
responsible not only for the currents in the skin layer of the ocean, but also for most of 
the circulation in the oceanic troposphere. 

We have long considered that the winds are responsible only for the current motion in 
the surface layer of about 200 meters. This depth is nothing but the “depth of the frictional 
influence” defined by Ekman. To explain the circulation in deeper parts of the troposphere, 
we had to assume a very strong convection current and slope current. Still we had a dis¬ 
tinct difference in the circulation patterns between the troposphere and stratosphere. 
These circumstances have made several problems very complicated. Defant (1928) defined 
the troposphere as the part of the ocean in which we can expect strong currents due to 
violent turbulence and convection. Still we can have violent convection in the seas of 
higher latitudes beyond the polar fronts which are no longer defined as troposphere. The 
conclusion that the drift currents penetrate into much deeper layers than D z is much in 
favor of the definition of the troposphere as the upper layer of the ocean in which strong 
currents are present. 

The explanation of the result that we can have a strong‘motion even in a layer several 
hundred meters deep might seem to be possible by assuming slope currents which would 
be produced as the effect of purely wind-driven water masses piled up against the land 
barriers. As a matter of fact, Ekman’s theory assumes no boundaries and a constant 
latitude. We can prove the existence of slope current in an ocean having boundaries partly 
or completely enclosing it. The slope current is uniform from the surface down to the 
bottom. This fact seems in favor of the theoretical result we have obtained. Still, we must 
give attention to the fact that the velocity of slope current is always inversely proportional 
to the depth of the sea. When the depth is large, as we see in the actual oceans, the slope 
current will not be strong enough to account for those large velocities we have obtained 
at the depth twice or three times as large as D z . 

We do not know an appropriate explanation of the theoretical result that the effect of 
winds can be felt at a depth several times as large as Ekman’s depth of frictional influence. 
It is to be hoped that someone may be able to solve this question satisfactorily in the near 
future. 

SUMMARY 

(1) A theory of the general circulation of water in the Pacific Ocean produced by the 
semipermanent wind system prevailing over this ocean is propounded. 

(2) The velocity is used to express the water motion which has formerly been explained 
by several authors in terms of mass transport. 

(3) The Pacific Ocean is considered to be a rectangular ocean extending from 60° S to 
60° N latitudes and from 0° to 120° longitude, and a zonal distribution of the wind system 
determined from actual observations has been assumed. 

(4) The effects of horizontal turbulence and the meridional variation of the Coriolis 
force have been taken into account. 

(5) The patterns of horizontal circulation are obtained in terms of streamlines for the 
sea surface and several deeper layers specified by the ratio z/D z where z is the geometrical 
depth below the surface and D z the depth of the frictional influence, a measure of vertical 
turbulence. 
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(6) Surface circulation has a pattern similar to that actually observed and does not 
differ much from Munk’s result obtained in terms of mass transport. We have very strong 
western currents and boundary vortices. 

(7) The magnitude of the Kuroshio and other western currents was computed from the 
distribution of the streamlines in each level. The velocity is inversely proportional to D z 
so that we can determine it by assuming an appropriate value for D z , A value of D z be¬ 
tween 50 m. and 150 m. gives values most reasonable and consistent with the actual ob¬ 
servations. 

(8) ' Subsurface circulations also show similar patterns except for a general decrease in 
motion as we go down into deeper layers. Still, it is remarkable that we have far stronger 
currents than expected from Ekman’s classical theory even at a depth much-larger than 
Ekman’s depth of frictional influence at which we can scarcely expect any motion. This 
seems to show us that the winds are responsible for most of the tropospheric motion of 
water. 
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TABLE 6 

Four Roots of the Equation: a 4 — 1082323+a- 2 — 14997110a- + 51883042+ = 0 (m = 2) 
For Different Values of tj. Refer to Equation (73) 


V 

a 

0 

T 

8 

0.1 

261.2089 

\ - 130.6217 

0.034595 

( -130.6217 



( - 200.88201 


( +200.88201 

0.2 

304.2429 

i - 152.1906 

0.138326 

\ -152.1906 



( - 161.71271 


( +161.71271 

0.3 

371.1853 

i - 185.7480 

0.310731 

\ -185.7480 



( - 7734621 


<+ 77.34621 

0.4 

454.0497 

- 363.0414 

0.55003 

- 91.55829 

0.5 

545.9024 

- 489.8041 

0.85179 

- 56.95000 

0.6 

642.5958 

- 603.9481 

1.20755 

- 39.85528 

0.7 

741.9602 

- 713.6338 

1.60417 

- 29.93054 

0.8 

842.8718 

- 821.2038 

2.02475 

- 23.69276 

0.9 

944.7270 

- 927.6127 

2.45104 

- 19.56531 

1.0 

1047.1851 

- 1033.3250 

2.86653 

- 16.72663 

1.1 

1150.0451 

- 1138.5916 

3.25872 

- 14.71229 

1.2 

1253.1819 

- 1243.5582 

3.61994 

- 13.24361 

1.3 

1356.5153 

- 1348.3155 

3.94678 

- 12.14657 

1.4 

1459.9923 

— 1452.9222 

4.23898 

- 11.30907 

1.5 

1563.5762 

- 1557.4175 

4.49829 

- 10.65703 

1.6 

1667.2415 

- 1661.8286 

4.72746 

- 10.14036 

1.7 

1770.9698 

- 1766.1750 

4.92960 

- 9.72442 

1.8 

1874.7475 

- 1870.4707 

5.10800 

- 9.38487 

1.9 

1978.5647 

- 1974.7263 

5.26555 

- 9.10440 

2.0 

2082.4135 

- 2078.9493 

5.40496 

- 8.86915 

2.5 

2601.9675 

- 2599.7504 

5.90334 

- 8.12046 

3.0 

3121.8036, 

- 3120.2640 

6.19651 

- 7.73611 

3.5 

3641.7746 

— 3640.6434 

6.38114 

- 7.51234 

4.0 

4161.8168 

- 4160.9508 

6.50416 

- 7.37016 

4.5 

4681.9004 

- 4681.2161 

6.58996 

- 7.27426 

5.0 

5202.0097 

- 5201.4554 

6.65206 

- 7.20632 

5.5 

5722.1359 

- 5721.6778 

6.69840 

- 7.15650 

6.0 

6242.2734 

- 6241.8885 

6.73386 

- 7.11876 

6.5 

6762.4187 

- 6762.0907 

6.76160 

- 7.08960 

7.0 

7282.5705 

- 7282.2877 

6.78369 

- 7.06649 

7.5 

7802.7263 

- 7802.4800 

6.80157 

- 7.04787 

8.0 

8322.8854 

- 8322.6689 

6.81623 

- 7.03274 

8.5 

8843.0470 

- 8842.8553 

6.82841 

- 7.02011 

9.0 

9363.2108 

- 9363.0397 

6.83861 

- 7.00973 

9.5 

9883.376 

- 9883.222 

6.84727 

- 7.00089 

10.0 

10403.542 

-10403.404 

6.85474 

- 6.99330 








Pacific Ocean Circulation — HlDAKA 


213 


TABLE 7 

Four Roots of the Equation: <r 4 — 1082323+<r 2 — 14997110«r + 103741765+ = 0 (m = 3) 
For Different Values of t?„ Refer to Equation (73) 


V 

a 

P 

y 

8 

0.1 

261.19836 < 

i - 130.63377 

! - 200.881841 

0.06917 j 

) -130.63377 
\ +200.88184 i 

0.2 

304.21377 l 

) - 152.24513 

! - 161.735111 

0.27648 j 

\ -152.24513 
! +161.735111 

0.3 

371.14543 < 

| - 185.88275 

! - 77.576941 

0.62007 < 

\ -185.88275 
!+ 77.576941 

0.4 

454.0086 

- 362.9384 

1.09300 

- 92.16317 

0.5 

545.8643 

- 489.7452 

1.67853 

- 57.79759 

0.6 

642.5616 

- 603,9042 

2.34715 

- 41.00454 

0.7 

741.9296 

- 713.5981 

3.05872 

- 31.39024 

0.8 

842.8444 

- 821.1731 

3.77053 

- 25.44186 

0.9 

944.7023 

- 927.5862 

4.44709 

- 21.56323 

1.0 

1047.1631 

- 1033.3018 

5.06562 

- 18.92693 

1.1 

1150.0246 

- 1138.5702 

5.61601 

- 17.07041 

1.2 

1253.1630 

- 1243.5387 

6.09746 

- 15.72172 

1.3 

1356.4978 

- 1348.2976 

6.51453 

- 14.71474 

1.4 

1459.9760 

- 1452.9056 

6.87418 

- 13.94458 

1.5 

1563.5610 

- 1557.4020 

7.18401 

- 13.34300 

1.6 

1667.2272 

- 1661.8141 

7.45125 

12.86434 

1.7 

1770.9563 

- 1766.1614 

7.68234 

- 12.47727 

1.8 

1874.7348 

- 1870.4579 

7.88287 

- 12.15980 

1.9 

1978.5526 

- 1974.7141 

8.05757 

- 11.89611 

2.0 

2082.4020 

- 2078.9377 

8.21039 

- 11.67465 

2.5 

2601.9583 

- 2599.7412 

8.74444 

- 10.96153 

3.0 

3121.7968 

- 3120.2572 

9.05081 

- 10.59044 

3.5 

3641.7687 

- 3640.6376 

9.24114 

- 10.37229 

4.0 

4161.8116 

- 4160.9456 

9.36693 

- 10.23297 

4.5 

4681.8958 

- 4681.2115 

9.45421 

- 10.13849 

5.0 

5202.0056 

- 5201.4513 

9.51716 

- 10.07142 

5.5 

5722.1322 

- 5721.6741 

9.56401 

- 10.02208 

6.0 

6242.2700 

- 6241.8851 

9.59981 

- 9.98471 

6.5 

6762.4151 

- 6762.6871 

9.62776 

- 9.95572 

7.0 

7282.5676 

- 7282.2848 

9.64999 

- 9.93278 

7.5 

7802.7236 

- 7802.4773 

9.66797 

- 9.91430 

8.0 

8322.8828 

- 8322.6663 

9.68271 

- 9.89921 

8.5 

8843.0446 

- 8842.8528 

9.69494 

- 9-88672 

9.0 

9363.2083 

- 9363.0372 

9.70520 

- 9.87627 

9.5 

9883.3736 

- 9883.2201 

9.71389 

- 9.86743 

10.0 

10403.5401 

-10403.4015 

9.72132 

- 9.85989 
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TABLE 8 

Four Roots of the Equation: o- 4 - 1082323+a 2 - 14997110a + 207804544+ = o (m = 4) 
For Different Values of r?. Refer to Equation (73) 


V 

a 

P 

y 

8 

0.1 

261.17780 < 

l - 130.658175 

1 - 200.8864001 

0.138550 < 

1 -130.658175 
! +200.8864001 

0.2 

304.15535 < 

1 - 152.354359 

1 - 161.7848571 

0.553368 < 

1 -152.354359 
! +161.7848571 

0.3 

371.06547 < 

1 - 186.151299 

| - 78.0307861 

1.237127 j 

1 -186.151299 
\ + 78.030786! 

0.4 

453.9263 

- 362.7315 

2.162989 

- 93.35781 

0.5 

545.7878 

- 489.6268 

3.271038 

- 59.43202 

0.6 

642.4929 

- 603.8161 

4.469333 

- 43.14513 

0.7 

741.8683 

- 713.5264 

5.657710 

- 33.99958 

0.8 

842.7895 

- 821.1121 

6.758457 

- 28.43583 

0.9 

944.6527 

- 927.5328 

7.73047 

- 24.85037 

1.0 

1047.1177 

- 1033.2539 

8.56385 

- 22.42764 

1.1 

1149.9834 

- 1138.5273 

9.26719 

- 20.72327 

1.2 

1253.1250 

- 1243.4993 

9.85663 

- 19.48229 

1.3 

1356.4627 

- 1348.2616 

10.35076 

- 18.55185 

1.4 

1459.9433 

- 1452.8723 

10.76555 

- 17.83659 

1.5 

1563.5304 

- 1557.3709 

11.11539 

- 17.27486 

1.6 

1667.1985 

- 1661.7851 

11.41206 

- 16.82546 

1.7 

1770.9293 

- 1766.1340 

11.66511 

- 16.46041 

1.8 

1874.7092 

- 1870.4321 

11.88225 

- 16.15935 

1.9 

1978.5283 

- 1974.6896 

12.06966 

- 15.90836 

2.0 

2082.3790 

- 2078.9146 

12.23233 

- 15.69674 

2.5 

2601.9398 

- 2599.7227 

12.79226 

- 15.00941 

3.0 

3121.7831 

- 3120.2435 

13.10803 

- 14.64763 

3.5 

3641.7570 

- 3640.6258 

13.30241 

- 14.43361 

4.0 

4161.8014 

- 4160.9354 

13.43018 

- 14.29618 

4.5 

4681.8867 

- 4681.2024 

13.51850 

- 14.20280 

5.0 

5201.9974 

- 5201.4431 

13.58205 

- 14.13631 

5.5 

5722.1247 

- 5721.6666 

13.62926 

- 14.08736 

6.0 

6242.2632 

- 6241.8783 

13.66528 

- 14.05018 

6.5 

6762.4092 

- 6762.0812 

13.69337 

- 14.02137 

7.0 

7282.5617 

- 7282.2789 

13,71571 

- 13.99851 

7.5 

7802.7181 

- 7802.4718 

13.73376 

- 13.98006 

8.0 

8322.8777 

- 8322.6612 

13.74855 

- 13.96505 

8.5 

8843.0398 

- 8842.8480 

13.76081 

- 13.95261 

9.0 

9363.2031 

- 9363.0320 

13.77110 

- 13.94218 

9.5 

9883.369 

- 9883.215 

13.77982 

- 13.93334 

10.0 

10403.535 

-10403.397 

13.78725 

- 13.92583 
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TABLE 9 

Four Roots of the Equation: + - 1082323*?V 2 - l4997110cr + 288890008+ = 0 {m = 5) 
For Different Values of rj . Refer to Equation (73) 


V 

a 

0 

y 

8 

0.1 

261.16177 J 

i - 130.677185 

! - 200.890771 

0.19260 J 

\ -130.677185 
• +200.890771 

0,2 

304.10979 < 

| - 152.43931 

| - 161.821521 

0.76882 < 

1 -152,43931 
\ +161.821521 

03 

371.00312 j 

1 - 186.35885 

1 - 78.38164! 

1.71458 < 

1 -186.35885 
!+ 78.381641 

0 A 

453.8619 

- 362.5695 

2.97958 

- 94.27196 

0.5 

545.7281 

- 489.5344 

4.45733 

- 60.65104 

0,6 

642.4394 

- 603.7475 

5.99960 

- 44.69154 

0.7 

741.8205 

- 713.4706 

7.46727 

- 35.81721 

0.8 

842.7467 

- 821.0646 

8.77347 

- 30.45556 

0.9 

944.6141 

- 927.4913 

9.88811 

- 27.01095 

1.0 

1047.0828 

- 1033.2171 

10.81805 

- 24.68375 

1,1 

1149.9513 

- 1138.4939 

11.58649 

- 23.04389 

1.2 

1253.0955 

- 1243.4691 

12.22046 

- 21.84683 

1.3 

1356.4353 

- 1348.2335 

12,74497 

- 20.94673 

1.4 

1459-9178 

- 1452.8463 

13.18123 

- 20.25277 

1.5 

1563.5066 

- 1557.3468 

13.54645 

- 19.70629 

1.6 

1667.1763 

- 1661.7625 

13.85429 

- 19.26805 

1.7 

1770.9083 

- 1766.1128 

14.11561 

- 18.91107 

1.8 

1874.6893 

- 1870.4120 

14.33896 

- 18.61630 

1.9 

1978.5095 

- 1974.6706 

14.53110 

- 18.36997 

2.0 

2082.3610 

— 2078.8964 

14.69740 

- 18.16196 

2.5 

2601.9254 

- 2599.7082 

15.26688 

- 17.48409 

3.0 

3121.7694 

- 3120.2297 

15.58610 

- 17.12578 

3.5 

3641.7452 

- 3640.6140 

15.78195 

- 16.91313 

4.0 

4161.7911 

- 4160.9251 

15.91042 

- 16.77647 

4.5 

4681.8775 

- 4681.1932 

15.99912 

- 16.68342 

5.0 

5201.9892 

- 5201.4349 

16.06288 

- 16.61715 

5.5 

5722.1172 

- 5721.6591 

16.11022 

- 16.56832 

6.0 

6242.2563 

- 6241.8714 

16.14639 

- 16.53126 

6.5 

6762.4032 

- 6762.0752 

16.17447 

- 16.50244 

7.0 

7282.5558 

- 7282.2730 

16.19684 

— 16.47963 

7.5 

7802.7126 

- 7802.4663 

16.21491 

- 16.46125 

8.0 

8322,8725 

- 8322.6560 

16.22972 

- 16.44623 

8.5 

8843.0349 

- 8842.8431 

16.24200 

- 16.43379 

9.0 

9363.1991 

- 9363.0280 

16.25230 

- 16.42337 

9.5 

9883.3649 

- 9883.2114 

16.26102 

- 16.41455 

10.0 

10403.5319 

-10403.3933 

16.26847 

- 16.40704 
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TABLE 10 

Four Roots of the Equation: o- 4 — 1082323+<r 2 — 14997110o- + 468733228+ = 0 0 = 6) 
For Different Values of rj. Refer to Equation (73) 


V 

a 

(3 

7 


8 

0.1 

261.1262 

\ - 130.71934 

0.312479 

5- 

130.71934 

0.2 

304.0087 

/ - 200.8931H 

\ - 152.62721 

1.245716 

( +200.8931 11 
\ -152.62721 

0.3 

370.8647 

l - 161.905761 
\ - 186.81402 

2.763347 

( +161.905761 
\ -186.81402 

0.4 

453.7192 

(- 79.148541 

- 362.2092 

4.74125 

} + 

79.148541 
96.25127 

0.5 

545.5957 

- 489.3292 

6.94392 

— 

63.21045 

0.6 

642.3205 

- 603.5949 

9.10051 

— 

47.82611 

0.7 

741.7145 

- 713.3466 

11.02080 

— 

39.38867 

0.8 

842.6516 

- 820.9590 

12.63321 

— 

34.32580 

0.9 

944.5284 

- 927.3991 

13.94741 

— 

31.07676 

1.0 

1047.0047 

- 1033.1347 

15.00642 

— 

28.87640 

1.1 

1149.8801 

- 1138.4198 

15.85925 

— 

27.31956 

1.2 

1253.0299 

- 1243.4015 

16.54939 

— 

26.17782 

1.3 

1356.3746 

- 1348.1714 

17.11211 

— 

25.31536 

1.4 

1459.8613 

- 1452.7887 

17.57493 

— 

24.64758 

1.5 

1563.4537 

- 1557.2930 

17.95898 

— 

24.11967 

1.6 

1667.1265 

- 1661.7121 

18.28046 

— 

23.69487 

1.7 

1770.8615 

- 1766.0655 

18.55182 

— 

23.34779 

1.8 

1874.6451 

- 1870.3674 

18.78268 

— 

23.06042 

1.9 

1978.4675 

- 1974.6283 

18.98052 

— 

22.81973 

2.0 

2082.3211 

- 2078.8563 

19.15124 

— 

22.61604 

2.5 

3.0 

3121.7487 

— 3120.2090 

20.05542 

_ 

21.59512 

3.5 

4.0 

4161.7756 

- 4160.9061 

20.38232 

_ 

21.25182 

4.5 

5.0 

5201.9767 

- 5201.4224 

20.53547 

_ 

21.08974 

5.5 

6.0 

6242.2460 

- 6241.8611 

20.61914 

_ 

21.00404 

6.5 

7.0 

7282.5470 

- 7282.2642 

20.66976 

_ 

20.95256 

7.5 

8.0 

8322.8648 

- 8322.6483 

20.70268 

_ 

20.91919 

8.5 

9.0 

9363.1903 

- 9363.0192 

20.72530 

_ 

20.89636 

9.5 

10.0 

10403.524 

-10403.385 

20.74146 

- 

20.88003 







